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摘 要 : 目前 土壤 盐 渍 化 是 全 球 重 要 的 环境 问题 , 探 明 于 田 绿洲 土壤 盐 渍 化 的 时 空 变化 规律 ,挖掘 雷达 遥感 探测 


土壤 盐分 的 优势 ,对 干旱 区 绿洲 的 土壤 盐 渍 化 时 空 变 化 进行 监测 评估 。 以 于 田 绿 洲 为 研究 区 ,基于 PALSAR-2、 
Sentinel-1 极 化 合成 孔径 雷达 数据 和 Landsat 8 OLI 等 多 源 数 据 集 , 租 选 雷 达 影像 的 最 优 后 向 散射 特征 与 主 成 分 分 析 
后 的 光学 影像 组 合 , 最 后 利用 随机 森林 方法 进行 图 像 分 类 ,定量 提取 于 田 绿洲 土壤 盐 渍 化 信息 ,对 土壤 盐 渍 化 时 空 
变化 进行 分 析 。 结 果 表明 :(1) 在 同时 使 用 随机 森林 分 类 方法 下 ,各 年 的 光学 影像 总 体 精度 平均 为 80.36% ,Kappa 


系数 平均 为 0.77; 光 学 影像 结合 雷达 影像 的 分 类 精度 比 光 


学 影像 分 类 精度 高 ,总体 精度 平均 为 85.62%,Kappa 系数 


平均 为 0.82,(2) 2015—2021 年 于 田 绿洲 产生 土壤 盐 渍 化 的 区 域 主要 分 布 于 研究 区 北部 的 绿洲 边缘 和 荒漠 交错 
带 。(3) 2015—2021 年 盐 渍 地 面积 年 均 变 化 量 为 -1120.55 hm2.a-, 变 化 率 为 -10.67%。 于 田 绿 洲 盐 渍 化 程度 总 体 呈 


下 降 趋势 , 盐 渍 化 以 轻 中 度 盐 渍 地 为 主 。 


关键 词 : 土壤 盐 涡 化 ; 合成 孔径 雷达 ; Landsat 8 OLI; 随机 森林 分 类 ; 时 空 变化 


土壤 盐 渍 化 是 全 世界 目前 共同 面临 的 生态 环 
境 问 题 , 盐 渍 化 问题 会 间接 损害 生态 系统 ,严重 阻 
得 干旱 、 半 干旱 区 当地 的 农业 生产 活动 和 经 济 的 可 
持续 发 展 "。 我 国 盐 溃 土壤 分 布 广 泛 且 各 种 盐 大 十 
所 占 面积 约 9.9x107 hma, 新 疆 的 盐 渍 化 土壤 总 面积 
约 8.5x10’ hm, 是 由 当地 气候 条 件 和 独特 的 地 形 地 
貌 所 造成 3。 利用 遥感 手段 对 绿洲 地 区 土壤 盐 渍 
化 时 空 分 布 和 土壤 盐 渍 化 程度 进行 监测 ,实现 当地 
自然 生态 环境 稳定 和 确保 绿洲 农业 生产 的 健康 发 
展 起 到 重要 作用 “1， 

遥感 技术 监测 区 域 范围 内 的 土壤 盐 渍 化 始 于 
20 世 纪 70 年 代 ”。 国 内 外 涉及 遥感 的 应 用 中 ,光学 
图 像 的 应 用 最 为 广泛 ,但 由 于 光学 遥感 属于 被 动 遥 
感 ,在 获取 数据 过 程 中 因 其 光谱 特性 易 受到 入 射 光 


二 进行 监测 ,但 日 间 成 像 和 受 云雨 影响 的 情况 下 ， 
盐 涡 化 信息 的 提取 依赖 光谱 特征 具有 局 限 性 ""。 
雷达 遥感 作为 主动 过 感 ,不 受气 候 和 剧 夜 影 
啊 , 具 备 全 天 时 、 全 天 候 的 成 像 特点 , 且 能 提供 光学 
遥感 不 包含 的 信息 。 融 合 了 微波 遥感 等 多 源 和 遥感 
数据 的 监测 方法 可 以 弥补 光学 遥感 数据 的 局 限 性 ， 
为 目标 检测 提供 光谱 和 纹理 特征 的 独特 组 合 ”。 目 
前 基于 微波 遥感 数据 监测 土壤 盐 渍 化 的 相关 研究 
中 ,学 者 们 通过 分 析 岁 像 后 向 散射 强度 .土壤 盐分 
和 介 电 常数 的 关系 ,发现 微波 影像 可 以 作为 土壤 盐 
分 监测 的 有 效 手段 "。 利 用 雷达 全 极 化 数据 ,分 析 
雷达 后 向 散 冉 系数 与 土壤 盐分 之 间 的 关系 ,证 明 结 
合 极 化 组 合 特征 的 微波 影像 能 对 不 同 程 度 的 盐 滇 
化 进行 区 分 。 使 用 光学 和 雷达 影像 多 源 数据 , 基 


源 和 天 气 条 件 影响 ,导致 土壤 信息 的 监测 分 类 精度 
降低 ,光学 遥感 利用 地 物 的 光谱 特征 对 盐 渍 化 土 
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化 监测 方面 的 可 行 性 ”。 

综 上 ,集成 多 源 光 学 与 雷达 数据 应 用 在 盐 渍 化 
监测 的 研究 已 有 一 定 成 果 ,但 对 盐 渍 化 时 空 变化 监 
测 的 相关 研究 较 少 。 对 此 ,本 文 探讨 利用 多 源 遥 感 
数据 和 野外 实测 数据 ,通过 提取 有 效 的 影像 特征 信 
息 进行 最 优 特征 组 合 , 对 人 研究 区 使 用 随机 森林 分 类 
方法 ,识别 土壤 盐 渍 化 的 分 布 并 分 析 结 果 , 人 研究 于 
田 绿 洲 土 壤 盐 渍 化 的 时 空 变化 特征 ,探索 微波 与 光 
学 数据 监测 干旱 区 绿洲 土壤 盐 碱 化 的 能 力 。 


1 研究 区 及 数据 来 源 


1.1 研究 区 概况 

F FA AR (36°30'~37°05'N , 81°09’ ~82°03’E) {iz 
Fie reg PAA) BLU ee. 
东 接 民 丰 县 , 西 邻 策 勒 县 , 南 与 西藏 自治 区 日 土 县 
相 接 ,绿洲 面积 约 3.95x10' km"), HAE A AdE 
低 , 垂 直 地 带 差异 明显 ,相对 高 差 4142 m, HI E 
腿 状 所 。 绿 洲 内 季节 性 河流 发 源 于 昆仑 出 的 克 里 
雅 河 属于 以 冰川 融雪 为 补给 的 内 陆 河 , 流 经 于 田 绿 
洲 , 消 失 在 塔克拉玛干 沙漠 腹地 ,当地 水 资源 缺乏 ， 
水 系 相对 封闭 中 。 该 区 域 属于 暧 温带 内 陆 干旱 启 
潢 气候 ,昼夜 温差 大 ,多 年 平均 气温 为 12.2 % ,最 高 
温度 39.4% ,最 低温 度 -17.6 % ;蒸发 作用 强 , 年 降水 
量 约 33.5 mm; 绿 洲 植被 种 类 贫乏 ,生物 群落 结构 单 
一 ,生态 系统 较 脆弱 ,土壤 盐 渍 化 现象 严重 |。 
1.2 数据 资料 与 处 理 
1.2.1 ERAH 雷达 数据 首先 选取 日 本 航天 局 发 
射 的 对 地 观测 卫星 (Advanced Land Observing Satel- 


81°20'E 81°30'E 
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37°0'N 
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lite,ALOS-2) 所 获取 成 像 时 间 为 2015 年 4 月 23 日 
的 (Phased Array type L-band Synthetic Aperture Ra- 
dar,PALSAR ) 全 极 化 数据 。 其 次 选择 了 欧洲 空间 局 
(European Space Agency, ESA ) 开 发 的 Sentinel-1 E 
星 , 成 像 日 期 为 2018 年 4 月 23 日 .2021 年 4 月 18 日 
的 Sentinel-1A SAR 双 极 化 数据 。 使 用 ENVI 5.3 $K 
件 中 的 SARscape 5.2.1 模块 对 影像 进行 多 视 处 理 、 
3 mx3 m 窗 口 精细 Lee 滤 波 、 地 理 编码 辐射 定 标 ,并 
重 采 样 为 20 m 分 辨 京 ,为 与 Sentinel-1 极 化 方式 保 
持 一 致 ,选择 VV、VH 极 化 方式 进行 极 化 组 合 。 

光学 遥感 数据 采用 Landsat 8 陆地 成 像 仪 (Oper- 
ational Land Imager, OLI) 数据 ,成像 时 间 分 别 为 
201544 H 26 H 2018 年 4 月 18 日 和 2021 年 4 月 10 
日 ,具体 参数 信息 见 表 1。 使 用 ENVI 5.3 软件 对 时 
序 Landsat 8 数据 进行 辐射 定 标 、 大 气 校正 、 影 像 配 
准 、 重 采样 为 20 m 分 辨 率 等 预 处 理 流程 。 
1.2.2 野外 数据 ”本 研究 根据 于 田 绿 洲 的 土壤 盐 污 
化 和 交错 带 分 布 空间 规律 ,于 2015 年 4 月 22 日 至 5 
月 7 日 和 2021 年 7 月 9 日 至 7 月 15 日 进行 野外 样品 
的 采集 ,筛选 出 74 个 有 代表 性 的 采样 点 ,分布 如 图 
la 所 示 ,每 个 采样 点 取 0~20 cem KE EIE .拍照 并 记 
录 土 壤 类 型 .周围 植被 类 型 建立 采样 点 图 片 库 便 
于 对 研究 结果 进行 精度 评价 。 将 样本 进行 自然 风 
干 . 人 研磨 和 使 用 1.5 mm 角子 过 得 ,使 土壤 样品 和 蒸 
饮水 混合 ,溶液 配制 为 1:5 的 土 水 比 ,使 用 电导 仪 在 
25 人 C 的 室温 下 测量 电导 率 。 通 过 构建 电导 率 和 总 
可 游 性 盐 之 间 的 方程 式 对 土壤 的 含 盐 量 进行 计算 
All pH {ERWE 
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yÈ: Bla F RX Band 7,G 4 Band 5,B J Band 2; 图 b 中 R 为 Band VH ,G X Band VV,B 为 Band VLH+VV。 下 同 。 
图 1 研究 区 影像 及 采样 点 


Fig. 1 Images and sampling points of study area 
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表 1 遥感 数据 参数 


Tab. 1 Parameters of remote sensing data 


卫星 传感器 获取 日 期 分 辩 率 /m 产品 类 型 极 化 方式 /波段 
PALSAR-2 2015-04-23 10x10 Level 1.1 VV. VH 
Sentinel- 1A 2018-04-23 10x10 GRD VV. VH 
Sentinel- 1A 2021-04-18 10x10 GRD VV.VH 
Landsat8 OLI 2015-04-26 30x30 LevellT 1~7 
Landsat8 OLI 2018-04-18 30x30 Level 1T 1~7 
Landsat8 OLI 2021-04-10 30x30 


Level lT 1~7 


2 研究 方法 


首先 获取 于 田 绿 洲 研 究 区 的 野外 实测 数据 、 光 
学 和 雷达 多 源 数据 。 然 后 将 光学 和 雷达 数据 以 多 
源 融 合 的 方式 ,通过 特征 提取 、J-M 距 离 分 析 和 随机 
和 森林 分 类 方法 ,对 人 研究 区 土壤 盐 污 化 进行 分 析 ,最 
后 对 研究 区 土壤 盐 渍 化 动态 变化 进行 监测 。 有 具体 
技术 路 线 如 图 2 所 示 。 
2.1 特征 提取 及 分 析 

遥感 影像 依赖 不 同 的 地 物 光 谱 信 息 区 分 地 物 类 
别 ,单一 的 特征 信息 无 法 提供 足够 的 区 分 信息 ”。 
将 不 同 的 特征 信息 进行 融合 ,融合 后 的 特征 数据 在 
一 定 程度 解决 了 单一 特征 信息 片面 问题 ,将 各 种 特 
征 描述 的 信息 聚合 ,可 提升 高 分 辩 率 影像 的 分 类 精 
度 ”。 

首先 将 预 处 理 后 的 光学 影像 进行 主 成 分 分 析 
(Principal Component Analysis, PCA) ,再 对 雷达 影像 
进行 后 向 散射 特征 提取 并 组 合 , 最 后 使 用 Pearson 相 
关 分 析 , 计 算 电 导 率 EC 值 与 后 向 散射 组 合 特征 间 的 
显著 性 。 根 据 图 3 筛选 出 最 优 后 向 散射 特征 组 合 。 
以 主 成 分 分 析 后 的 光学 影像 为 基础 ,将 各 分 类 特征 
进行 琶 加 组 合 ,建立 分 类 特征 数据 集 ( 表 2) ,用 同样 
的 方法 得 到 不 同年 份 下 的 相同 分 类 特征 数据 集 。 
2.2 J-M 距 离 

为 定量 研究 不 同 分 类 特征 下 地 物 间 的 可 分 离 
程度 ,对 训练 样本 使 用 J-M 上 距离 评估 地 物 间 的 可 分 
离 性 ,J-M 上 距离 的 数值 范围 在 0~2 之 间 , 若 两 种 地 物 
的 距离 值 越 趋 近 2, 说 明 两 种 类 型 的 可 分 离 性 越 
高 ”。J-M 计 算 方法 如 下 : 

J-M=2(1-e”) 

式 中 :8B 表示 巴 氏 距离 。 


于 田 绿洲 研究 区 


' r 请 “ALOS—2 PALSAR-2, 


多 视 处 理 
EAM 处 十。 外表 


可 分 离 性 评价 
随机 森林 分 类 


于 多 源 数据 的 盐 || 土壤 盐分 空间 
化 识别 精度 分 析 | | 动态 变化 分 析 


土壤 盐分 时 间 
序列 变化 分 析 


a ey et a i it ee el 


y 
土壤 盐 渍 化 动态 变化 监测 


Fig.2 Workflow of study 


Sie 


计算 结果 如 表 3, 方 式 三 光学 主 成 分 特征 和 
SAR 影像 的 极 化 组 合 特征 铸 加 下 ,重度 盐 添 地 与 其 
他 地 物 的 J-M 距 离 最 佳 ,表明 各 地 物 的 可 分 离 性 在 
光学 和 雷达 数据 的 组 合 下 可 区 分 性 较 好 。 

2.3 随机 森林 分 类 

将 随机 森林 (CRF) 算 法 用 于 土壤 盐 渍 化 的 分 
类 。 随 机 森林 以 集成 分 类 顺 的 形式 进行 分 类 ,分 类 
顺 由 众多 决策 树 组 合 而 成 ,具有 训练 速度 快 、 泛 化 
性 能 强 和 模型 稳定 等 优点 ”。 随 机 森林 方法 对 于 
集成 多 源 和 遥感 数据 ,具有 多 个 分 类 特征 的 数据 集 可 
以 有 效 提高 分 类 精度 “。 徐 乔 等 ”对 极 化 SAR 影 
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注 :EC 为 土壤 电导 率 ;* 在 0.05 级 别 ( 双 尾 ) 相 关 性 显著 。 
图 3 Pearson 相关 性 和 矩阵 


Fig. 3 Pearson correlation matrix diagram 


像 进行 分 类 ,利用 随机 森林 分 类 算法 实验 得 到 高 精 
BE .光滑 上 且 分 类 层次 连续 的 分 类 结果 。 本 文 结合 文 
献 [27-30] 中 的 方法 ,根据 特征 数 来 对 决策 树 的 个 
数 进行 调试 。 


3 结果 与 分 析 


3.1 基于 多 源 数据 的 盐 渍 化 识别 精度 分 析 
为 对 比 不 同 分 类 特征 条 件 下 地 物 的 识别 效果 ， 
将 各 年 的 PCA 光谱 特征 影像 和 PCA+ 后 向 散射 特征 


影像 进行 分 类 结果 对 比 。 因 方式 二 ,雷达 影像 进行 
多 次 分 类 后 ,总体 精 度 平均 为 67.07% ,难以 达到 区 
分 地 物 的 目的 , 故 排除 方式 二 。 根 据 乔 木 等 ”对 新 
疆土 壤 盐 渍 化 程度 分 级 标准 ,将 地 物 分 为 裸 地 、 植 
被 .水体 、. 轻 中 度 盐 渍 地 和 重度 盐 涡 地 5 类 ( 轻 中 度 
盐 渍 地 为 有 盐 生 植 被 覆盖 ,pH 值 在 7.1~9.5; 重 度 盐 
渍 地 为 无 植被 覆盖 pH 值 在 9.5 以 上 )。 本 研究 以 高 
分 辩 率 影像 ( 奥 维 互动 地 图 ) 和 野外 图 片 库 作为 目 
视 解 译 的 参考 ,选取 训练 样本 和 验证 样本 。 对 不 
同 的 训练 样本 均 使 用 相同 的 验证 样本 进行 分 类 后 
的 评价 ,便于 科学 对 比分 析 分 类 后 的 验证 结 

根据 表 4 中 分 类 精度 可 知 ,方式 一 的 总 体 精度 
平均 为 80.36% ,Kappa 系数 平均 为 0.77。 方 式 三 的 
总 体 精 度 平均 为 85.62% , Kappa 系数 平均 为 0.82。 
结合 表 3 和 表 4, 随 机 森林 分 类 的 精度 说 明 多 源 数 据 
特征 组 合 的 分 类 效果 是 要 优 于 光学 遥感 数据 。 雷 
达 影 像 的 后 向 散射 系数 对 土壤 盐 涡 化 的 介 电 常 数 比 
较 敏感 ,土壤 含 盐 量 直 接 影响 介 电 常数 的 变化 ”, 利 
用 这 一 特性 将 雷达 和 光学 影像 结合 ,可 以 弥补 雷达 
影像 噪声 对 地 物 识 别 的 干扰 ,增强 后 向 散射 系数 对 
盐 渍 地 的 识别 从 而 提高 分 类 精度 。 

为 详细 对 比 不 同 数据 源 分 类 的 结果 ,以 分 类 精 
度 较 好 的 2021 年 为 例 , 选 择 于 田 绿 洲 具 有 代表 性 的 
[工区 域 进行 局 部 特征 的 对 比 。 如 图 4 所 示 , 区 域 
1 、 工 是 于 田 绿 洲 两 个 “牛角 边缘 ,沙漠 与 盐 澳 地 


表 2 分 类 特征 数据 集 组 合 


Tab. 2 Classification feature data set combination 


序列 组 合 方式 简称 波段 备注 
I PCA 光谱 特征 方式 一 3 - 
Il 后 向 散射 特征 方式 二 4 VV 
VV+VH 
VV2+HVH 
(VV2+VHVH 
亚 PCA+ 后 向 散射 特征 FRE 7 一 
RI 典型 地 物 的 J-M 距离 
Tab.3 J-M distance of typical ground objects 
样本 J-M 测 度 
地 物 数据 集 一 一 一 一 一 
轻 中 度 盐 涡 地 PRHE 植被 水 体 
重度 盐 渍 地 方式 一 1.710 1.928 1.762 1.924 
方式 二 1.274 1.249 1.456 1.684 
方式 三 1.893 1.899 1.973 2.000 
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表 4 分 类 精度 验证 


Tab.4 Classification accuracy verification 


oe 2015 4E 2018 年 2021 年 
Z| x 
总 体 精度 /96 Kappa 系数 总 体 精度 /% Kappa 系数 总 体 精度 /9% Kappa 系数 
方式 一 79.55 0.75 80.32 0.77 81.21 0.79 
方式 三 85.19 0.81 84.59 0.80 87.07 0.84 
81°20'E 81°30'E 81°40'E 
Z z I 
人 F 
cS œ 
on on 
I 
Z Z 
© © 
D D 
$ A mm 植被 m 水 体 加 裸 地 


81°20'E 


81°30'E 81°40'E 


注 :图 a 和 图 c 为 PCA 光 谱 特 征 分 类 结果 ;图 b 和 图 d 为 PCA+ 后 向 散射 特征 分 类 结果 。 
图 4 随机 森林 分 类 结果 对 比 


mm epee mm 重度 盐 渍 化 


Fig. 4 Comparison chart of random forest classification results 


ZO De Sah, CE eh oP , PEP BE A BBE Eb Tot 
地 交错 带 十 分 明显 。 区 域 1 体现 了 绿洲 水 体 和 盐 
涡 地 的 过 渡 交 错 分 布 ,PCA 光谱 特征 数据 的 分 类 结 
果 中 (图 4a) ,两 侧 盐 渍 地 存在 错误 分 类 的 散 点 状 植 
被 区 域 ,PCA+ 后 向 散射 特征 的 分 类 结果 较为 有 效 的 
减少 了 零碎 斑点 的 出 现 , 体 现 盐 涡 化 交错 带 的 边界 
特征 (图 4p)。 在 区 域 开 地 物 由 裸 地 沙漠 .重度 盐 渍 
地 和 轻 中 度 盐 活 地 逐渐 过 渡 到 水 体 区 域 ,图 4c 和 图 
4d 区 域 在 植被 和 中 轻 度 盐 省 地 的 区 分 上 存在 略微 
的 不 同 ,因为 植被 比例 在 15%~70% ,光学 影像 中 地 
表 反 射 率 可 能 表现 为 植被 和 其 他 地 物 的 混合 ”。 根 
据 野 外 采样 发 现 耐 盐 植 被 如 芦苇 、 树 柳 和 盐 生 草木 
可 以 在 轻 中 度 盐 渍 化 土壤 中 生长 ,这 些 耐 盐 植被 会 
间接 影响 光学 影像 中 盐 渍 化 土壤 的 光谱 反射 率 ”。 
PCA+ 后 向 散射 特征 影像 具有 的 后 向 散射 系数 可 以 
透 过 植被 冠 层 ,对 地 物 的 纹理 信息 进行 捕捉, 将 其 
结合 光学 遥感 的 光谱 信息 ,在 一 定 程度 上 增强 影像 
的 分 离 度 ,消除 植被 的 干扰 ,将 轻 度 盐 渍 地 的 稀 牙 
植被 进行 区 分 ,降低 其 对 地 物 分 类 的 影响 中 。 获 取 
植被 覆盖 下 的 盐 渍 土 含 盐 量 信息 ,还 需 依靠 微波 散 
射 模型 来 反 演 土壤 含 盐 量 ”。 因 此 ,将 被 动 的 光学 


遥感 和 主动 的 微波 遥感 数据 综合 是 目前 监测 土壤 
盐 渍 化 较为 理想 的 方法 ,未 来 研究 中 有 必要 结合 土 
坏 湿 度 .地 表 粗 糙 度 和 土壤 介 电 常数 来 更 好 的 解释 
土壤 盐分 的 分 布 情况 。 
3.2 土壤 盐分 空间 动态 变化 分 析 

从 分 类 结果 可 知 ( 图 5) ,于 田 绿洲 盐 渍 化 主要 
发 生 在 北部 西南 部 ,绿洲 边缘 和 荒漠 交错 带 区 域 ， 
这 与 相关 学 者 的 研究 一 致 ”。 盐 渍 化 区 域 中 , 重 
度 盐 渍 化 土壤 主要 在 研究 区 北部 的 两 个 “牛角 ”之 
间 分 布 ,以 及 克 里 雅 河流 域 的 下 游 区 域 , 呈 重度 盐 
渍 化 向 塔克拉玛干 沙漠 的 过 渡 形 式 。 在 绿洲 近郊 
轻 中 度 盐 涡 化 土壤 分 布 在 植被 和 重度 盐 渍 化 土壤 
之 间 ,呈现 条 带 状 分 布 ,主要 原因 是 绿洲 内 排 碱 渠 
中 的 水 未 经 处 理 排放 到 了 绿洲 边缘 ,导致 绿洲 边缘 
轻 中 度 盐 涡 化 程度 较 严 重 ”。 重 度 盐 渍 地 与 轻 中 
度 盐 渍 地 分 布 区 域 交 错 分 布 , 流 经 的 区 域 因 地 下 水 
位 通常 较 高 ,伴随 地 表 强 烈 的 蒸发 作用 ,土壤 盐 渍 
化 积累 作用 显著 ”。 

对 人 研究 区 表层 土壤 盐 渍 化 空间 动态 变化 分 析 ， 
提取 各 年 同一 时 期 盐 渍 化 的 面积 变化 ( 表 5) ,于 田 
绿洲 2015 年 的 盐 渍 化 面积 占 总 面积 的 40.12% ,其 
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81°20'E 81°30'E 81°40'E 
(a) 2015 年 


81°20'E 81°30'E 81°40'E 
(b) 201844 


81°20'E 81°30'E 81°40'E 
(c) 2021 年 
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mm 植被 m 水 体 国 裸 地 国 轻 中 度 盐 渍 化 m 重 度 盐 渍 化 
图 5 各 年 份 图 像 分 类 结果 


Fig. 5 Classification results of images by year 


表 5 201S 一 2021 年 地 物 分 类 面积 变化 
Tab.5 The area change of classified features from 2015 to 2021 


、 2015 年 2018 年 2021 年 
分 类 类 别 
面积 /hm 比例 /9% 面积 /hm? 比例 /9% 面积 /hm? 比例 /% 
水 体 4839.01 3.09 7652.32 4.88 5553.04 3.52 
植被 61061.10 38.94 68645.04 43.78 71506.76 45.31 
裸 地 27970.34 17.84 25340.48 16.16 24533.92 15.54 
轻 中 度 盐 渍 地 41518.64 26.48 36462.24 23.25 32359.16 20.51 
重度 盐 涡 地 21411.11 13.65 18700.08 11.93 23847.28 15.12 
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中 , 轻 中 度 盐 渍 化 区 域 占 26.48% ,重度 盐 渍 化 区 域 
占 13.65% 。2018 年 盐 渍 化 面积 占 总 体面 积 的 
35.18% , 其 中 轻 中 度 和 重度 盐 涡 化 区 域 各 占 
23.25% 11.93%。2021 年 盐 渍 化 区 域 占 总 体面 积 的 
35.63% , 其 中 轻 中 度 和 重度 盐 涡 化 区 域 各 占 
20.51% 、15.12%。2015 一 2021 年 于 田 绿 洲 的 土壤 盐 
渍 化 程度 总 体 呈 逐年 波动 下 降 趋 势 , 盐 渍 化 程度 主 
要 以 轻 中 度 盐 渍 地 为 主 。 轻 中 度 盐 渍 地 的 面积 变 
化 从 2015 一 2021 年 呈 持 续 下 降 趋 势 ,面积 占 比 逐渐 
减少 ,植被 面积 呈 上 升 趋势 ;重度 盐 涡 地 的 面积 呈 
先 减少 后 增加 的 趋势 ;绿洲 表层 水 体 变 化 趋势 不 明 
显 , 但 土壤 水 分 是 土壤 盐分 运动 和 变化 的 主要 驱动 
因素 ,间接 影响 当地 土壤 盐分 的 空间 格局 变化 “1。 
3.3 土壤 盐分 时 间 序 列 变化 分 析 

为 反映 于 田 绿 洲 盐 渍 土 在 2015 一 2021 年 的 迁 
移 变 化 ,进行 到 置 分 析 得 到 盐 渍 地 时 空 变 化 分 布 
(图 6) ,由 于 PALSAR-2 影 像 范围 的 限制 ,2015 年 的 
研究 仅 在 光学 影像 的 分 类 基础 上 进行 。 使 用 桑 基 
图 直观 可 视 化 的 体现 了 于 田 绿 洲 各 地 类 之 间 的 结 
构 特 征 和 转 入 转 出 关系 (图 7)。 结 果 表 明 2015 一 
2021 年 于 田 绿洲 的 主要 5 个 乡镇 区 域 的 耕地 周边 ， 
大 部 分 区 域 处 于 盐 渍 化 程度 降低 的 变化 趋势 。 但 
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6 2015—2021 年 盐 渍 地 变化 分 布 
Fig. 6 Distribution of saline change from 2015 to 2021 


TENE AR GES AR PA HB AS BO ET EHE FEE Ek Tae 1 i 
的 现象 , 因 绿洲 内 部 土壤 在 排 盐 的 改良 措施 下 盐分 
含量 较 低 , 当 灌溉 和 排放 失衡 时 ,土壤 表层 下 的 可 
溶性 盐 离 子 随 毛管 水 运 移 到 地 表 富 积 , 造 成 耕地 形 
成 次 生 盐 渍 化“*”。 盐 溃 化 程度 增加 的 希 吾 勒 乡 和 
斯 也 元 乡 北 部 区 域 ,根据 实地 考察 和 DEM 影像 ,发 
现 区 域 属 于 盆地 地 形 ,地 下 水 的 水 位 和 矿 化 度 相对 
较 高 。 北 部 海拔 高 度 相 对 于 南部 低 约 200 m, H. 
于 燥 的 气候 导致 水 分 蒸发 ,加 重 了 土壤 盐 渍 化 的 
形成 。 


202303.00108v1 


chinaXiv 


1# 肖 森 天 等 :基于 光学 和 雷达 多 源 


图 7 2015—2021 年 地 物 类 型 转移 网 
Fig. 7 Feature type transfer map from 2015 to 2021 


由 研究 区 土壤 盐分 时 间 序 列 变化 统计 结 

OP yeep a et 
-6723.31 hm2.a ,变化 率 为 -10.68% ,说 明 盐 渍 化 程 
度 总 体 处 于 降低 趋势 。 其 中 轻 中 度 盐 渍 化 年 均 变 
化 量 为 -1526.58 hm2.a ,变化 率 为 -22.06% , 呈 降 低 
的 变化 趋势 ;重度 盐 渍 化 的 年 均 变 化 量 为 406.03 
hm2.a ,变化 率 为 11.38% ,重度 盐 溃 化 程度 的 变化 
过 程 是 较为 缓慢 的 增加 趋势 ,说 明 重度 盐 涡 化 土壤 
的 改良 较为 困难 ,这 与 前 人 对 于 田 绿 洲 不 同 程度 盐 
渍 化 土壤 变化 的 预测 结果 基本 一 致 *。 


4 结论 


以 于 田 绿 洲 为 研究 对 象 , 采 用 野外 实测 数据 、 

学 数据 和 雷达 数据 结合 的 方式 ,通过 特征 提取 和 
te al 
林 分 类 。 得 到 以 下 主要 结论 : 

(1) 对 比 2015 年 .2018 年 和 2021 年 分 类 结果 
EE ie tachi aan ncaa 
类 的 精度 ,得 到 光学 影像 的 PCA 光谱 特征 的 总 体 精 
度 平均 为 80.36% ,Kappa 系数 平均 为 0.77。PCA+ 后 
向 散射 特征 的 组 合 总 体 精度 平均 为 85.62% ,Kappa 
系数 平均 为 0.82。 精 度 评价 结果 表明 使 用 随机 森林 
分 类 方法 在 该 分 类 特征 下 分 类 精度 较 高 。 说 明 在 
光学 影像 分 类 的 基础 上 加 入 雷达 影像 是 未 来 盐 渍 
化 监测 中 较为 有 效 的 方法 手段 。 

(2) 通过 对 于 田 绿 洲 土壤 盐分 空间 动态 变化 分 
析 ,得 出 2015 一 2021 年 于 田 绿洲 产生 土壤 盐 渍 化 的 
区 域 主 要 分 布 于 研究 区 北部 的 绿洲 边缘 和 荒漠 交 
错 带 。 盐 涡 地 的 面积 变化 为 轻 中 度 盐 渍 地 呈 逐 渐 


遥感 
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降低 趋势 ,于 田 绿 洲 的 植被 面积 逐年 增加 ,重度 盐 
渍 地 为 先 减少 后 增加 的 波动 趋势 ,于 田 绿 洲 内 的 河 
流水 体面 积 变化 不 显著 。 

(3) 对 于 田 绿 洲 土 壤 盐 涡 化 时 序 变化 分 析 , 数 
据 表明 2015—2021 年 盐 涡 化 程度 总 体 呈 下 降 趋 势 ， 
盐 涡 地 面积 年 均 变 化 为 -1120.55 hm a, ZAER 
为 -10.68% ,说 明 近 7 a 土壤 盐 渍 化 现象 得 到 较为 有 
效 控制 和 改良 。 
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Spatial and temporal analysis of soil salinity in Yutian Oasis by combined 


optical and radar multi-source remote sensing 


XIAO Sentian'*, Ilyas NURMEMET"’, Nuerbiye MUHETAER”, 
Zhao Jing'’, Adilai ABULAITI’* 


(1. College of Geographical and Remote Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 
2. Key Laboratory of Oasis Ecology, Xinjiang University, Urumqi 830046, Xinjiang, China) 


Abstract: Soil salinization is currently a prominent global environmental problem. Spatiotemporal variation of 
soil salinization in Yutian Oasis was explored, and the advantages of radar remote sensing in excavated soil 
salinity were investigated to monitor and evaluate the temporal and spatial variabilities of soil salinization in arid 
oasis. Based on Phased Array type L-band Synthetic Aperture Radar 2 (PALSAR-2), Sentinel- 1 polarimetric 
synthetic aperture radar data, and Landsat 8 Operational Land Imager (Landsat 8 OLI) multi-source dataset, the 
optimal backscattering characteristics of radar images were selected, and optical images were combined based on 
principal component analysis. Finally, the random forest method was used to classify the images. Quantitative 
extraction of soil salinization information in Yutian Oasis and the spatial and temporal variation of soil 
salinization were analyzed. Results showed that (1) under the random forest classification approach, the total 
accuracy of optical images in each year was 80.36%, and the kappa coefficient was 0.77. The classification 
accuracy of optical images combined with radar images was higher than that of optical images, the total accuracy 
was 85.62%, and the kappa coefficient was 0.82. (2) From 2015 to 2021, the area of salinized soil in Yutian Oasis 
was mainly distributed in the north of the study area, interlaced edge of the oasis, and desert. (3) The average 
annual variation of saline land area from 2015 to 2021 was -1120.55 hm’ - a`', and the change rate was — 10.67%. 
The salinization of Yutian Oasis generally showed a downward trend and was mainly mild to moderate saline 
land. 

Keywords: soil salinization; synthetic aperture radar; Landsat 8 OLI; random forest classification; spatial and 


temporal variation 


